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SECTION  I 


INTRODUCTION 

There  is  a  continuing  interest  on  the  part  of  the  U.  S.  Army  in 
fin-stabilized  vehicles  designed  to  fly  at  low  altitudes  and  hypersonic 
velocities.  Various  missiles  operate  in  this  flight  regime  as  do  the 
high  energy,  sabot-launched  projectiles  (KE  penetrators)  used  to  penetrate 
heavy  armored  vehicles. 

The  flight  environments  (surface  shear,  heat  transfer  and  pressure) 
needed  to  support  design  activities  traditionally  have  been  obtained  by 
correlating  laboratory  or  range  data  and  extrapolating  to  flight  conditions. 
It  is  now  recognized  that  finite  difference  methods  may  supplement  or  sup¬ 
plant  the  traditional  empirical  approaches  and  offer  both  a  better  under¬ 
standing  of  the  flowfields  and  a  potentially  significant  reduction  in 
design  costs.  The  use  of  finite  difference  methods  for  this  application 
is  being  assessed  and  explored  via  a  number  of  ongoing  studies  at  various 
research  organizations. 

Most  current  flowfield  models  consist  of  one  of  several  alternative 
sets  of  flowfield  equations,  a  perfect  gas  equation  of  state  and  simple 
correlations  or  tables  for  the  properties  (viscosity,  specific  heat  and 
Prandtl  Number).  Such  perfect  or  ideal  gas  formulations  are  valid  for 
air  flowfields  provided  that  (1)  the  Mach  numbers  are  not  high  enough 
to  cause  significant  dissociation,  (2)  there  are  no  condensibles  in  the 
flow  and  (3)  there  is  no  source  of  non-air  elements  within  the  flow.  Real 
gas  effects  can  be  important  and  should  be  included  in  the  formulation 
when  any  of  these  constraints  are  not  satisfied. 

The  real  gas  terms  associated  with  dissociation  enter  the  flowfield 
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mode]  via  the  equation  of  state  and  the  thermodynamic  and  transport 
properties.  These  terms  were  formulated  and  included  in  an  automated 
solution  procedure  for  the  parabolized  Navier-Stokes  equations  on  a 
previous  study  for  AFWAL.1  This  solution  procedure  is  used  in  the  present 
study  to  assess  the  significance  of  real  gas  effects  on  the  predicted 
flowfields  about  typical  fin-stabilized  vehicles  flying  at  hypersonic 
velocities  and  low  altitudes. 

Flowfield  predictions  have  been  obtained  at  (or  near)  standard  at¬ 
mospheric  conditions  for  three  Mach  Numbers  (M-6,  8  and  10).  The  vehicle 
considered  consists  of  a  spherically  tipped  small  half  angle  cone  joined 
to  a  cylindrical  afterbody.  Predictions  were  made  both  with  and  without 
real  gas  effects  to  allow  a  direct  assessment  of  real  gas  effects  on 
the  aerodynamic  and  aerothermodynamic  environments. 

The  following  sections  discuss  (1)  the  real  gas  parabolized  Navier- 
Stokes  flowfield  model,  (2)  the  calculations  performed  as  part  of  the 
present  study  and  (3)  the  conclusions  reached  as  a  result  of  the  present 
study. 
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SECTION  II 
REAL  GAS  TERMS 


Real  gas  terms  enter  the  parabolized  Navier-Stokes  relations2  via 
the  equation  of  state,  the  shock  fitting  logic,  the  S  matrix  and  the 
Jacobian  matrices  A,  B,  C,  M,  and  Ay,  where  the  terminology  has  been 
defined  previously.2  Perfect  gas  approximations  are  also  currently  em¬ 
ployed  in  the  calculation  of  the  speed  of  sound.  Modifications  to  the 
present  formulation  for  each  of  these  areas  will  be  discussed  in  the 
following  paragraphs. 

1.  Equation  of  State 

The  conservation  equations  are  written  in  terms  of  an  energy  variable 

.  ★ 
eT  defined  as  a  total  (internal  plus  kinetic)  energy  per  unit  volume. 

In  terms  of  commonly  employed  thermodynamic  variables,  eT  can  be  written 

as 

eT  -  pe  ♦  p[u2  +  v2  +  w2]/2  (1) 

where  e  is  the  specific  internal  energy  (energy  per  unit  mass).  For  a 
perfect  gas  pe  =  p/(y-l),  and  Eq.  (1)  reduces  to  the  form  usually  seen 
in  the  literature 

Cj  *  P  i  +  p[u2  +  v2  +  w2]/2  (2) 


★ 

Most  authors  designate  the  total  internal  per  unit  volume  as  e.  The 
change  to  e^  here  was  done  to  avoid  confusion  with  the  specific  internal 
energy,  needed  to  describe  real  gas  flows. 


For  a  real  gas  having  constant  elemental  composition. 


e  «  *(P,p)  £  (3) 

p 

where  <j> ( p ,p )  is  a  function  which  must  be  determined  numerically.  Sub¬ 
stituting  Eq.  (3)  into  Eq.  (1)  yields 

eT  *  0(p,p)  P  +  p[u2  +  v2  +  w2]/2  (4) 

which  now  properly  includes  real  gas  effects.  Note  that  Eq.  (4;  is 
nonlinear  in  pressure  in  contrast  to  Eq.  (2).  Methods  for  obtaining 
pressure  from  Eq.  (4)  will  be  discussed  subsequently. 

The  partial  derivatives  of  pressure  with  respect  to  the  elements 
(q)  of  the  U  matrix  (i.e.,  p,  pu,  pv,  pw,  e^.)  will  be  needed  in  forming 
the  Jacobian  matrices.  It  is  convenient  to  obtain  them  now  by  taking  the 
total  derivative  of  Eq.  (4).  Thus, 

deT  ■  (*  *  "  Ip)  dp  -  |g2  +  ,2  *  "2  -  57  (5) 

♦  ud(ou)  +  vd(pv)  +  wd(pw) 

The  partial  derivatives  of  pressure  are  easily  extracted  from  Eq.  (5)  as 


!?  -  *°-42  *  *2  *  "2  - 2p  IfM*  *  p  $)  (6a> 

-f£-j  .  -u/(*  ♦  p  f|)  (6b) 

sfouT  "  ',/(d  *  p  Ip)  (6c) 

iuT)  ■  -W(‘  *  P  Ip)  (6d) 
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which  completes  the  formulation  of  the  new  pressure  partial  derivatives 
for  the  Jacobian  matrices. 

2.  The  $  Matrix 

The  elements  of  the  S  matrix  represent  the  terms  on  the  right  side 
of  the  conservation  equations.  The  perfect  gas  formulation  is  correct  for 
the  continuity  and  momentum  equations  for  any  thin-layer  flow.  However, 
the  energy  equation  includes  a  flux  divergence  term  on  the  right  side 
which  includes  two  key  relations  applicable  only  to  perfect  or  ideal 
gases.  Revision  of  the  formulation  starts  with  a  generally  valid  form  of 
the  flux  divergence  term  (FDT)  which  can  be  written  as 


where  Cartesian  coordinates  are  employed  and  q  is  the  flux,  y  is  the 
coordinate  in  the  wall  normal  direction  and  pr,  Up,  and  L  are  reference 
values  of  density,  velocity  and  length  used  to  nondimensionalize  the 
equation.  In  the  perfect  gas  formulation,  the  relations 


(8) 


(9) 


are  used  to  reduce  Eq.  (7)  to  the  form 


(10) 
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where  nondimensional  variables  are  used  exclusively  and  are  defined  as 
follows: 


®T  -  eT /orUr 
t  =  t/(L/Ur) 

K  =  k/k„ 


Re  =  prUrL/wr  (Reynolds  number) 


Pr  =  ur  Cp  / kr  (Prandtl  number) 
ur  =  a. 

a  =  a/Ur  (sound  speed) 


where  C  is  taken  as  the  specific  heat  at  a  constant  pressure  for  the 
pf 

perfect  or  ideal  gas.  This  variable  will  be  redefined  subsequently  fo 
real  gas  as  the  frozen  specific  heat  (hence  the  f  subscript). 

The  real  gas  analog  of  Eq.  (8)  must  include  chemical  changes  in  the 
gas  and  the  associated  effects  of  diffusion.  The  heat  flux  becomes 


Q 


I 


species 


jihi 


(ID 


where  the  summation  is  taken  over  all  the  species  present  in  the  gas,  j. 

i.L 

is  the  mass  flux  of  the  i  species  and  h^  is  the  species  enthalpy 
(energy  per  unit  mass)  and  is  well  approximated  as  being  a  function  of 
temperature  only  =  h^(T). 

The  mass  flux  terms  appearing  in  Eq.  (11)  can  be  formulated  by 
introducing  Fick's  law  and  employing  the  binary  diffusion  approximation. 


where  V  is  the  binary  diffusion  coefficient  and  is  the  mass  fraction 
of  the  i*^  species.  Substituting  Eq.  (12)  and  j\  in  Eq.  (11)  yields  the 
heat  flux  q  as 


T  V  3k-i 

2  W 


3y  species 
which  is  further  simplified  to 


(13) 


3T 


q  .  -k  —  -  pD\—  -  C  — 


ay 


ah 


aT 


ay  Pf  ay 


(14) 


In  obtaining  Eq.  (14)  from  Eq.  (13),  the  following  definitions  and 
identities  have  been  employed: 


hiIkiN 

Cp  5  2  ki  cp.  (frozen  specific  heat) 


y  k3 

Z* .  ay 

species  J 


&  2  kicp- 

J  species  Ki 


^  ak .  ^  ^  ah. 

2  hiiF  =  ly  Z  kihi  -  2  kiir 


species 


species  species 


(15) 

(16) 


(17) 


(18) 


The  Schmidt  number  Sc  defined  as 


is  usually  employed  to  reduce  Eq.  (14)  to  the  form 


(19) 


(20) 
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which  is  in  a  form  suitable  for  nondimensional ization  and  subsequent 
evaluation. 


Turbulence  terms  can  be  included  in  the  heat  flux  formulation  by 
adding  turbulence  contributions  to  the  thermal  conductivity  and  diffusion 
coefficient.  Equation  (13)  with  turbulence  becomes 


q 


-(k  *»‘h  w[i£ 


(21) 


where  eH  is  the  eddy  conductivity,  is  the  eddy  diffusivity  and  both 
are  related  to  the  eddy  viscosity  eM  via  the  relations 


CH  =  Pr 

(22) 

_  €_H_ 

CD  =  Sct 


where  Pr^  and  Sct  are  the  turbulent  Prandtl  and  Schmidt  numbers,  respec¬ 
tively.  The  generalized  form  of  Eq.  (20)  which  includes  turbulence 
effects  is  readily  obtained  as 


This  equation  can  be  simplified  into  the  commonly  used  effective  conduc¬ 
tivity  form,  i .e. 

«  ■  "keff  !y 

for  flows  having  constant  elemental  composition  and  negligible  pressure 
gradients  in  the  y  direction.  Generally,  these  are  valid  constraints 
in  the  region  near  a  wall  where  heating  rate  is  important  provided  there 
is  no  source  of  foreign  elements  at  the  wall  (via  blowing,  transpiration, 
ablation,  etc.).  Thus,  the  enthalpy  gradient  can  be  written  as 
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Hi  -  c  — 

ay  "  pe 

where  C  is  the  equilibrium  specific  heat.  Substituting  this  relation 


into  Equation  (23)  yields  a  definition  for  the 

c  =  iS- 
(eff  Sc 


k.«  -  fc  !c.  +  (I7  -  1)  Cpf  *  5^  Cpe  +  (  Pr^  -  0  Cpf 


Lpe 


Sc 

W 


which  is  in  its  final  form  prior  to  the  introduction  of  nondimensional 
variables. 

The  additional  nondimensional  variables  needed  are  defined  as 


C  =  Cn  T  lut 

Pf  Pf  r  r 


Cp  =  Cp  T  /u2 
He  Ke  r  r 


T  =  T/T 

r 

eM  h  pr  cM/pr 
h  =  h/uf 


Substituting  these  relations  into  Eq.  (24)  and  rearranging  terms  yields 


fdt  -  h  *  %y v  [fe  (#  -  >) ♦  ?  -  ■)  1 ; Epf !  j  1251 

Finally,  Eq.  (25)  is  transformed  to  generalized  spatial  coordinates  to  yield 


FDT 


-*  *  d  *  d>  * 
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l^c  ^ct /  ^e 


‘)J  EpfJ  \ 


(26) 


to  be  used  in  place  of  the  existing  term,  i.e.. 


(24) 
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FDT  -  [uj;  +  z*  *  cfjo”1  K  Pr-1  (y  .  i)-l(a2)J^ 


where  the  Reynolds  number  has  been  factored  out  of  both  Eqs.  (26)  and  (27) 
and  appears  elsewhere  in  the  conservation  equation. 

Equation  (26)  has  two  terms  of  the  form  8^  d  6^  'K  each  of  which  is 


differenced  as 


^£+1  +  -  'J'fc)  -  +  -  'J'jui 

o..2 


which  is  the  same  difference  equation  presently  used  to  evaluate 
Eq.  (27).  The  linearity  of  Eq.  (28)  is  retained  by  lagging  the  property 
evaluation  —  i.e.,  using  properties  evaluated  at  the  previous  time  (or 
space)  step.  This  allows  the  procedure  developed  for  the  perfect  gas 
formulation  to  be  employed  here. 

3.  Jacobian  Matrices 


The  first  three  of  the  Jacobian  matrices  to  be  considered  are  of 


the  form 


A  =  J'1^  0(IU)/3U] 
B  .  J'1 Ca  e(V,n)/aU] 
6  «  J'1 [3  0{W,S)/8U] 


(29a) 


(29b) 


(29c) 


where 


pUa  ♦  B„P 


0  (a,  B) 


pVa  +  8  P 


pWa  +  B_P 


(eT+  p)o  -  BtP 
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Real  gas  effects  impact  the  partial  derivative  terms  only  via  terms  of  the 
form 


which  are  readily  obtained  from  the  relations  presented  in  Eq.  (6). 

A 

The  NR  matrix  is  the  Jacobian  of  the  real  gas  analog  of  the  S 
matrix  (which  we  will  relabel  as  SR).  The  formulation  is 

Mr  *  J-l  £aSR/aO  ]  (31) 


where 


The  top  four  elements  in  SR  [SR(1),  SR(2),  SR(3),  and  SR(4)]  are 


identical  to  the  top  four  elements  in  §  and  need  not  be  modified.  The 

A 

bottom  element  SR ( 5 )  contains  the  real  gas  contribution  to  the  flux 
divergence  and  is  significantly  different  from  the  corresponding  element 
S(5). 
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Derivatives  of  the  various  thermodynamic  and  transport  properties 
are  required  in  order  to  support  evaluation  of  the  a$R/aO  elements 

A 

appearing  in  the  MR  matrix.  The  properties  in  question  include  T,  b,  u. 
Sc,  C  ,  and  Pr.  The  turbulent  Prandtl  and  Schmidt  numbers  will  be 
taken  to  be  constants  (Prt  .  0.9,  Sct  -  0.9),  while  the  eddy  viscosity 
is  dealt  with  separately. 

Let  x(p,  p)  represent  any  one  of  the  properties  of  interest  (i.e. 
T,  h,  u,  •  •  .)•  The  chain  rule  is  used  to  obtain  the  total  derivative 


dx  -iii  dp  +  dp 
3p  '  p  3p'p 

allowing  dp  to  be  obtained,  i.e. 


"p- (dl  -f?iP  d»)/(l?ip) 

and  substituting  this  result  into  Eq.  (5)  to  obtain 

deT  =  Rdx  -  u2  +  v2  +  w2  -  2p  ~  +  2R  —1 
i  (_  3 p  3p|  2 


+  u  d(pu)  +  v  d(pv)  +  w  d(pw) 


where 


t  *  -|t 

3  in 


Eq.  (34)  can  be  used  to  generate  all  of  the  partial  derivatives  for  each 
of  the  properties.  Thus, 


ill 

3 p  'pU 


b  +  v2  +  w2  -  2p  —  I  +  2R  — 1”| 

,  pv,  pw,  eT  w  6p  ap  ip  3p  IpJ 


3pu  ip,  pv,  pw,  eT  r 


ax  .  _  _ 

apv'p,  pu,  pw,  e-j-  r 


(38) 


ax  ■  _  _  w 

apwlp,  pu,  pv,  ey  r 


(39) 


ax  |  _  1^ 

aej'p,  pu,  pv,  pw  r 


(40) 


which  completes  the  formulation  of  the  property  derivatives. 

The  partial  derivatives  aS^/aU  are  obtained  by  substituting 
finite  difference  relations  for  the  quantities  u ’  ,  v  ,  w  ,  K  ,  T  then 

^  ^  ^  S  S 

using  the  chain  rule  on  the  resulting  expression.  This  is  a 
straightforward  operation  yielding  a  lengthy  result  which  will  not  be 
presented  here. 

4.  Shock  Fitting 

Shock  fitting  procedures  are  employed  in  the  evaluation  of  both 
unsteady  and  steady-state  supersonic  flows  to  determine  the  position  of 
the  bow  shock  wave  and  some  of  the  flow  and  thermodynamic  variables  behind 
it.  Initially,  the  pressure  behind  the  bow  shock  (p2)  is  obtained  by 
explicitly  integrating  a  combination  of  the  equations  of  motion  from  the 
wall  to  the  shock.  Knowing  p2  allows  the  shock  velocity  (unsteady 
flows)  or  slope  (steady-state  flows)  to  be  determined  from  the 
Rankine-Hugoniot  relations.  The  variable  (shock  velocity  or  slope) 
obtained  from  this  evaluation  is  then  used  to  obtain  the  new  shock 
position  by  propagating  the  shock  in  the  c -direction  along  constant  £  and 
n  lines.  Real  gas  effects  enter  the  formulations  via  the  Rankine-Hugoniot 
relations.  The  details  follow  for  both  the  unsteady  and  steady-state 
shock  fitting  procedures. 
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In  the  unsteady  flow  formulation,  the  shock  is  allowed  to  move 
outward  away  from  the  body  with  a  velocity  q*s  in  the  shock  normal 
direction.  The  sign  conventions  for  qg  and  the  other  important  velocity 
components  are  shown  in  Fig.  1.  The  shock  velocity  is  given  by  the  relatic 


q  =  u  +  u, 

HS  SSj  1 


(41) 


where  u^  is  the  component  of  the  freestream  velocity  vector  q*i 
directed  toward  the  shock  along  the  shock  normal,  andTTss  is  the 
velocity  vector  the  shock  would  have  if  q^  were  set  to  zero  while 
holding  P2 / Pi  fixed.  The  vector  -u^$  can  also  be  viewed  as  the 
freestream  velocity  component  in  the  shock  normal  direction  required  to 
maintain  a  steady-state  shock  wave  of  pressure  ratio  p^/pj. 

The  unknown  velocity  vectors  are  obtained  by  simple  manipulation  of 
the  vector  relations  starting  with  Uj  wh i ch  is  obtained  by 


Uj  «  dj  • 


(42) 


where  ns  is  the  unit  vector  in  the  shock  normal  direction.  It  follows  that 
“l  -  ifVx  '  *  Vy  3  *  "l«l  ^  <43> 


where 


m  =  {x}  +  x}  +  x}) 
x  y  2 7 


1/2 


(44) 


and  Uj,  vj,  wj  are  the  components  of  q*!  and  cx,  cu,  Cz  are 

the  coordinate  metrics.  The  variable  uss-j  can  be  obtained  directly  for 

a  perfect  gas  from  the  relation 


USS,I  * 


(45) 
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and  its  known  direction.  The  replacement  formulation  including  real  gas 
effects  requires  two  steps.  First,  the  density  behind  the  shock  is 
obtained  from  the  relation 


[Pj  +  (1  +  Zi2)  p2 
P2  +  (1  +  29j)  Pj 


(46) 


which  is  the  real  gas  analog  of  the  perfect  gas  relation 


p2  =  Pj 


( Y  +  l)Po  +  (Y  -  1)  Pi 


(y  +  lJPi  +  (y  -  l)  r 


(47) 


which  also  appears  in  the  perfect  gas  formulation.  Having  obtained  both 
the  pressure  and  density  behind  the  shock,  |uss^  |  is  readily  obtained 
from  the  Rankine-Hugoniot  relation 


or  in  final  form 

■“ssj  ■-lu«s1l[cx  S  ^  kl( 


(«8) 


m 


(49) 


The  velocity  component  u$$  can  also  be  evaluated  behind  the  shock  (i.e., 

Tf  ).  This  variable  is  interpreted  as  the  shock  normal  component  of 
2 

velocity  required  to  maintain  a  steady-state  shock  of  pressure  ratio 
P2/Pj.  It  is  evaluated  from  the  relation 


ss. 


Po  ss 


1 


(50) 
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or 


P1|USS] 
ss^  ni 


''x**',5'  C2  k] 


a 


The  components  of  velocity  tangent  to  the  shock  wave  (v^,  are 
given  by 


V1  =  v2  *  ql  -  U1 


(52) 


or 


-  ui  (‘  -  r)5 *  vi(‘  -  *  “i^  -  ;r)' 


(53) 


which  completes  the  evaluation  of  the  components  needed  to  evaluate"^. 

The  velocity  components  are  added  in  the  manner  illustrated  in 
Fig.  1  to  yield 


q2  *  uss2  +  v2  +  qs 


(54) 


which  becomes 


f?  - - -  u  +  q*i  -  u,  +  "u*.  +  u ,  =  q*i  +  IT  (l - 

^  o  2  s  s  ^  X  X  ssj  X  X  sSj  y  p  ^  y 


(55) 


after  substituting  in  Eqs.  (41),  (49),  and  (51).  The  components  of  q2 
are  readily  evaluated  from  Eqs.(48)  and  (54)  as 


u2  =  Uj  +  B  cx/ m 

v2  =  V1  +  B  V"1 

w2  *  W1  +  B 


(56) 
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where 


.  h)  pnhS  A  -H) 

\Vpl/  \p2  •  pl/  \  p2/ 


(57) 


Eq.  (56)  is  identical  to  the  analogous  relation  in  the  perfect  gas 
formulation.  Equation  (57)  is  significantly  different  in  that  B  is 
related  to  perfect  gas  variables  as 


(58) 


Equations  (47),  (56),  and  (57)  represent  the  changes  to  the 
unsteady  flow  shock  fitting  procedure  associated  with  real  gas  effects. 

The  remaining  formulation  is  unchanged  from  the  perfect  gas  formulation  as 

3 

described  by  Kutler  and  Pedelty. 

For  the  PNS  application,  the  velocity  component  normal  to  the  shock 
is  obtained  from  a  real  gas  version  of  the  Rankine-Hugoniot  relation 


J 


(Px  -  P2) CP1  +  +  P2] 


(59) 


which  is  used  in  place  of  the  equivalent  perfect  gas  relation 


/( y  ♦  1)  P?  +  (y  -  1 

i  “.i  -  J - ^ — 


-  1)  p. 


(60) 


The  metrics  at  the  shock  are  obtained  either  from  coordinate 
transformation  relations  (cx,  r?x,  £x)  or  from  numerical 
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differentiation  (for  cyf  c2,  n  n2,  £yi  S2)  in  the  crossflow 
plane.  The  formulation  for  these  relations  is  unchanged  by  the 
introduction  of  real  gas  effects,  e.g. 

Wj  *  vlV  *  tMuV z  *  v/  -  (“l  -  Cv)] 


C.  V 

(Uj  -  Uj) 


n  =  C 
X  X 


S&lfii' 


A  typical  point  on  the  shock  is  propagated  using  the  following 


equations: 


where 


zj+1  -  ac| 

yJ+1  =  yJ  +  yJ 


ZK  *  (Vx  “  Vy)/J 
y?  =  (Vz  -  Vx,/J 


The  density  component  behind  the  shock  wave  is  calculated  from  the  real 


gas  relation 


@2*0 


fPj  +  (1  +  2^2)p2"| 
l[p?  ♦-(!■*  J*jTpjJ 


which  is  used  in  place  of  the  perfect  gas  relation 


(y  ♦  l)Po  +  (y  -  1)P, 


1  Pi  +  Y 


and,  finally,  the  velocity  components  behind  the  shock  are  found  from 


where 


u2  *  U1  +  Bss  V* 

vr'i*  Bss  Vm 

"2  =  W1  *  Bss  V 


Bss  =  l“l I (1  ■  °l/p2) 


which  completes  the  formulation. 

5.  Implementation  in  PNS  Code 

The  flowfield  procedures  obtain  real  gas  properties  from  tabulations 
stored  in  memory.  The  data  in  the  tables  is  generated  prior  to  and  sepa¬ 
rate  from  the  flowfield  calculation.  Tabular  data  for  air  was  obtained 

4  5 

from  well  established  compilations  ’  and  is  presently  included  in  the 
* 

code.  Data  for  environments  other  than  air  can  be  generated  using  a 

6 

procedure  discussed  elsewhere  .  The  material  presented  herein  discusses 
the  data  requirements,  the  format  of  the  tables,  the  searching  procedures 
employed,  and  the  interpolation  relations  used  to  obtain  data  between 
table  entries. 

The  independent  variables  employed  are  log  p  and  log  p/p  ,  where 

10  1 0  O 

pQ  is  the  density  of  air  at  standard  conditions  (pQ  =  1.225  x  10"3  gm/cm3). 
The  dependent  variables  include  four  thermodynamic  variables,  three  trans¬ 
port  properties,  and  two  partial  derivatives  of  transport  properties. 

The  thermodynamic  properties  appear  in  the  format  ♦  ,  ^T>  Cn  /C. 

P  P  Pf  P 

The  transport  properties  appear  in  the  format  y,  Pr,  and  Sc.  In  addition, 

*  The  real  gas  viscosity  data  never  reduces  to  the  Sutherland  formula 
exactly,  so  small  differences  will  always  exist  between  real  gas  and 
perfect  gas  predictions. 
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entropy  is  related  to  pressure  and  density  via  the  relation 


where  the  effective  specific  heat  ratio  y  is  also  tabulated  as  a 
function  of  log  and  log  — .  For  air  the  reference  condition 

10  P0  10  p0 

S 

was  set  to  23.5372. 

Each  of  the  tabulated  properties  is  evaluated  at  the  last  marching 
step  as  known  functions  of  pressure  and  density,  then  frozen  and  treated 
as  being  independent  of  x  (or  ?)  in  the  evaluation  of  the  equations  at 
the  current  marching  step.  The  approach  is  successful  because  each  of 
the  tabulated  properties  is  a  slowly  varying  function  of  axial  distance. 


SECTION  III 


CALCULATIONS 

Calculations  have  been  carried  out  for  a  single  model  geometry 
flying  at  three  interesting  flight  conditions.  The  results  obtained 
allow  an  assessment  to  be  made  of  the  effects  of  the  real  gas  equation 
of  state  and  transport  properties  on  the  prediction  of  the  aerodynamic 
and  aerothermodynamic  environments. 

1.  Model  Geometry 

The  model  consists  of  a  cone-cylinder  with  an  8  degree  half  angle 
conical  tip,  0.075-inch  radius  spherical  nose  bluntness,  1.4  inch  diam¬ 
eter  body  and  an  overall  length  of  11.20  inches  (eight  caliber). 

2.  Flight  Conditions 

The  flight  conditions  for  each  of  the  three  cases  considered  are 
given  in  Table  1. 


Table  1 

Flight  Conditions 


Mach  Number 

Angle  of  Attack  (deg.) 

Free  Stream 

Static  Temp.  (°R) 

Static  Pressure  (atm) 

Wall  Temperature  (°R) 

Reynolds  No.  (in_1xl0”6) 


Case  Number 


1 

2 

3 

6 

8 

10 

2 

2 

2 

530 

530 

530 

1 

1 

1 

530 

530 

530 

3.615 

4.82 

6.025 

The  surface  geometry  was  generated  using  the  FLAG=2  option  (see 
Ref.  2),  plus  a  surface  geometry  file  having  45  points  in  the  circum¬ 
ferential  half-plane  and,  finally,  used  in  conjunction  with  a  virtual 
orgin  situated  .45865  inch  upstream  of  the  nose  (used  to  generate  the 
spherical  cap  and  the  conical  surface).  All  output  from  the  computer 
is  referenced  to  the  virtual  origin,  not  the  nose  of  the  model. 

4.  Starting  Conditions 

The  blunt  body  code  of  Kutler  et  al3  was  used  to  generate  a  flowfield 
solution  over  the  spherical  cap  of  the  model  for  each  case  considered. 

The  output  from  each  solution  consists  of  two  planes  of  data  to  be  used 
to  start  PNS  code  calculations.  For  each  of  the  cases  considered,  the 
downstream  data  plane  was  situated  at  a  distance  x  =  .0587  inch  from 
the  nose  which  is  just  upstream  of  the  tangency  point  at  x  =  .06456  inch 
from  the  nose.  At  the  starting  plane,  the  surface  slope  is  at  an  angle 
of  12.55°  to  the  axis,  so  that  most  of  the  very  rapid  expansion  associated 
with  the  nose  cap  is  over. 

The  Kutler  et  al3  blunt  body  code  is  a  perfect  gas  code  which  was  used 
only  because  no  alternative  real  gas  code  of  this  type  was  available. 
Consequently,  the  starting  solutions  were  always  somewhat  inconsistent 
with  the  real  gas  PNS  solutions.  The  sensitivity  of  the  Case  3 
solution  to  this  inconsistency  was  assessed  as  part  of  the  calculation 
matrix. 

5.  Downstream  Solutions 


An  initial  set  of  solutions  was  calculated  to  establish  the  validity 
of  the  prediction  procedure.  Our  objective  was  to  assess  the  sensitivity 


of  the  solutions  to  the  arbitrary  inputs  which  must  be  adopted  in  making 
calculations.  The  following  issues  are  of  particular  concern  and  were 
addressed: 

(1)  Does  either  the  position  of  the  starting  plane  of  data  and/or 
the  interpolations  needed  to  start  the  calculations  have  any 
lasting  effect  on  the  nature  of  the  downstream  solution? 

(2)  It  is  necessary  to  introduce  arbitrary  transition  regions 
to  gradually  change  the  flowfield  from  perfect  gas  to  real 
gas  and/or  from  laminar  to  turbulent  flow.  Do  the  features 
of  the  transition  region  selected  have  any  lasting  effect 
on  the  nature  of  the  downstream  solution?  The  transition 
region  was  never  allowed  to  extend  beyond  .1787  inch  down¬ 
stream  from  the  nose. 

(3)  Damping  and/or  smoothing  parameters  must  be  selected  and  em¬ 
ployed  to  stabilize  calculations.  Do  they  have  any  lasting 
effect  on  the  nature  of  the  downstream  solutions? 

(4)  The  distance  from  the  wall  to  the  first  grid  point  off  the 
wall  (DS)  must  be  selected  arbitrarily  and  must  be  small 
enough  to  maintain  accuracy  but  large  enough  to  maintain 
stability.  Can  suitable  values  be  found?  An  initial  DS 

of  .0006  inch  was  selected  for  laminar  flows  while  .0004 
inch  was  selected  for  turbulent  flows. 

It  was  found  that  the  selection  of  the  damping  and  smoothing  parameters 
have,  by  far,  the  most  important  effect  on  the  nature  of  the  downstream 
solution.  Indeed,  the  selection  of  excessive  damping  parameters  (par¬ 
ticularly  and  eg,  Ref  2)  were  found  to  change  the  qualitative  features  of 


29 


solutions.  Values  of  =  .05  and  eg  =  .10  yielded  acceptable  solution'-, 
both  in  terms  of  accuracy  and  stability,  and  were  selected  for  the  final 
calculations.  Changing  the  other  arbitrary  inputs  tended  to  move  the 
solutions  around  a  little.  Fortunately,  such  changes  tended  to  have  the 
same  kinds  of  effects  on  both  the  perfect  gas  and  the  real  gas  solutions 
and,  consequently,  were  not  of  much  concern  for  the  purposes  of  the 
present  study. 

Flowfields  over  slender,  spherically  capped  cones  have  been  studied 
by  a  number  of  investigators.  The  results  of  Harris7  are  presented  in 
Figure  2  as  an  example  where  a  =  0,  R^  =  .5  inch,  =  8  and  the  conical 
half  angle  is  7  degrees.  As  shown,  the  pressure  drops  rapidly  arouid  the 
nose  cap,  overexpands,  then  recovers  slowly  to  the  conical  pressure.  The 
recovery  to  the  conical  pressure  can  be  measured  in  tens  of  nose  diameters, 
an  effect  associated  with  the  strong  entropy  layer  in  the  shock  layer. 
Increasing  (or  decreasing)  either  the  angle  of  attack  or  the  conical  half 
angle  decreases  (or  increases)  the  distance  required  to  recover  the  conical 
pressure. 

Axial  pressure  distributions  obtained  in  the  present  study  are  pre¬ 
sented  in  Figures  3,  4  and  8.  They  all  display  the  expected  qualitative 
features.  Moreover,  the  M^  =  8  pressure  distributions  given  in  Figure  8 
are  in  acceptable  agreement  with  the  Figure  2  results  considering  the 
differences  in  Reynold's  Number,  angle  of  attack  and  cone  half  angle. 

The  =  10  solutions  are  presented  in  Figures  3  to  6  in  the  forebody 
region.  It  was  felt  that  real  gas  effects  would  be  present  here  if  they 
are  going  to  be  present  anywhere.  An  examination  of  the  starting  plane 
of  data  showed  that  real  gas  effects  were  clearly  important  there.  Real 
gas  pressures  differed  by  as  much  a  25%  from  their  perfect  gas  counterparts. 


However,  by  the  time  the  flow  is  allowed  to  expand  and  cool  through  the 
transition  region  and  progress  downstream  far  enough  to  establish  a  valid 
PNS  solution,  the  effect  is  never  more  than  a  few  percent.  Pressure  dis¬ 
tributions  for  the  two  pitch  planes  are  given  for  laminar  flow  in 
Figure  3.  Pressure  distributions  for  the  two  pitch  planes  are  given  for 
turbulent  flow  in  Figure  4.  Similar  effects  are  shown  for  heat  transfer 
in  Figures  5  and  6. 

The  fact  that  the  heat  transfer  rates  and  pressure  distributions  are 
changed  very  little  does  not  imply  that  the  structure  of  the  shock  layer 
is  virtually  unchanged,  however.  To  illustrate  the  point,  the  temperature 
distribution  across  the  shock  layer  for  x  =  2.053  inches  at  the  wind¬ 
ward  pitch  plane  is  shown  in  Figure  7.  As  shown,  the  perfect  gas  temper¬ 
atures  are  about  450  °R  higher  than  the  real  gas  temperatures  at  the  maxi¬ 
mum  temperature  of  the  entropy  layer.  However,  both  perfect  gas  and  real 
gas  distributions  become  nearly  identical  at  both  shock  front  and  wall 
boundaries. 

To  complete  the  calculation  matrix,  the  M  =  6  and  M  =  8  cases  were 

OO  CO 

run  and  are  presented  in  Figures  8  and  9  along  with  the  entire  =  10 
results.  The  real  gas  effects  on  the  =  10  is  small,  while  the  effects 
of  real  gas  on  the  other  two  solutions  cannot  be  separated  from  uncertain¬ 
ties  in  the  solution  procedure  and  can  be  viewed  as  negligible. 


31 


SECTION  IV 


CONCLUDING  REMARKS 

On  the  basis  of  the  present  study,  the  following  conclusions  have 
been  reached  regarding  the  real  gas  effects  on  the  prediction  of  flow- 
fields  about  Army  projectiles: 

o  There  appears  to  be  a  moderately  significant  effect  on  the 
nose  of  the  vehicle  where  temperatures  and  pressures  are  high 
and  molecular  dissociation  is  enhanced.  However,  this  conclu¬ 
sion  is  tentative  since  the  PNS  code  used  in  the  present  study 
is  not  applicable  to  the  nose  cap  region  of  the  flow, 
o  Real  gas  effects  change  the  heating  rates  and  pressure  distri¬ 
butions  by  a  few  percent.  Even  this  small  effect  disappears 
after  1  or  2  inches  of  axial  length, 
o  Real  gas  effects  are  important  in  determining  the  interior 
structure  of  the  shock  layer  where  the  temperatures  can  be 
changed  significantly. 


(a)  Sign  Convention  for  Velocity  Vectors  (b)  Components  of  Parti  cl 
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Figure  1.  Shock  Wave  Notation 
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Figure  4.  Pressure  Distribution  for  Case  3  (M  =  10)  Showing  Real  Gas 
Effects  for  Turbulent  Flow. 
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Figure  9.  Real  Gas  Circumferential  Pressure  Distribution  at  y  r  11.2  Inches  for  Cases 
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